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Abstract-This paper examines the joint problem of detec- 
tion and identification of a sudden and unobservable change 
in the probabiUty distribution function (pdf) of a sequence 
of independent and identically distributed (i.i.d.) random 
variables to one of finitely many alternative pdf 's. The ob- 
jective is quick detection of the change and accurate infer- 
ence of the ensuing pdf. Following a Bayesian approach, a 
new sequential decision strategy for this problem is revealed 
and is proven optimal. Geometrical properties of this strat- 
egy are demonstrated via numerical examples. 

I. Introduction 

Consider a sequence of i.i.d. random variables Xi, X2, . . ., 
taking values in some measurable space {E,£). The common 
probability distribution of the X's is initially some known prob- 
ability measure Pq on {E, £), and then, at some unobservable 
disorder time 9, the common probability distribution changes 
suddenly to another probability measure for some unobserv- 
able index /i S A4 = {1,...,A/}. The objective is to detect 
the change as quickly as possible, and, at the same time, to iden- 
tify the new probability distribution as accurately as possible, so 
that the most suitable actions can be taken with the least delay. 

This problem can be viewed as the fusion of two fundamen- 
tal areas of sequential analysis: change detection and multi- 
hypothesis testing. In traditional change detection problems, 
there is only one change distribution. Pi; therefore, the focus 
is exclusively on detecting the change time. Whereas, in tradi- 
tional sequential multi-hypothesis testing problems, there is no 
change time to consider Instead, every observation has com- 
mon distribution P^ for some unknown /i, and the focus is ex- 
clusively on the inference of /i. Both of these subproblems have 
been studied extensively. For recent reviews of these areas, we 
refer the reader to [1] and [2] and the references therein. 

However, the joint problem involves key trade-off decisions 
not taken into account by separately applying techniques for 
these subproblems. While raising an alarm as soon as the 
change occurs is advantageous for the change detection task, 
it is undesirable for the identification task because waiting 
longer provides more observations for inferring the change dis- 
tribution. Likewise, the unknown change time complicates 
the identification task, and, as a result, adaptation of existing 
sequential multi-hypothesis testing algorithms is problematic. 
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Decision strategies for the joint problem have a wide array 
of applications, such as fault detection and isolation in indus- 
trial processes, target detection and identification in national de- 
fense, pattern recognition and machine learning, radar and sonar 
signal processing, seismology, speech and image processing, 
biomedical signal processing, finance, and insurance. However, 
the theory has not been broadly developed. Nikiforov [3] pro- 
vides the first results for this problem, showing asymptotic op- 
timality for a certain non-Bayesian approach, and Lai [4] gen- 
eralizes these results through the development of information- 
theoretic bounds and the application of likelihood methods. In 
this paper, we follow a Bayesian approach to reveal a new op- 
timal strategy for this problem and we describe an accurate nu- 
merical scheme for its implementation. 

In Sec. II we formulate precisely the problem in a Bayesian 
framework, and in Sec. Ill we show that it can be reduced to an 
optimal stopping of a Markov process whose state space is the 
standard probability simplex. In addition, we establish a simple 
recursive formula that captures the dynamics of the process and 
yields a sufficient statistic fit for online tracking. 

In Sec. IV we use optimal stopping theory to substantiate 
the optimality equation for the value function of the optimal 
stopping problem. Moreover, we prove that this value function 
is bounded, concave, and continuous on the standard probabil- 
ity simplex and that the optimal stopping region consists of M 
non-empty, convex, closed, and bounded subsets. Also, we con- 
sider a truncated version of the problem that allows at most N 
observations from the sequence of random measurements. We 
establish an explicit bound (inversely proportional to N) for the 
approximation error associated with this truncated problem. 

In Sec. V we show that the separate problems of change 
detection and sequential multi-hypothesis testing are solved as 
special cases of the overall joint solution. We illustrate some 
geometrical properties of the optimal method and demonstrate 
its implementation by numerical examples for the special cases 
AI ~ 2 and AI = 3. Specifically, we show instances in which 
the M convex subsets comprising the optimal stopping region 
are connected and instances in which they are not. Likewise, 
we show that the continuation region (i.e., the complement of 
the stopping region) need not be connected. We refer the reader 
to [5] for complete proofs of the results. 

II. Problem statement 

Let (n,!F,F) be a probability space hosting random vari- 
ables e : n ^ {0,1,...} and fi : ft M = {1,...,M} 
and a process X = (X„)„>i taking values in some measurable 



space {E, £). Suppose that for every t > l,i G Ai, n > 1, and 
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for some given probabiHty measures Pq, Pi, . . 
known constants po G [0, 1], p G (0, 1), and > 0, i G 
such that i^i + • • • + i^A/ = 1, where x A y ^ min{.T, y} and 
xM y = maxjx,?;}. Namely, 6' is independent of /.j; it has a 
zero-modified geometric distribution with parameters po and p 
in the terminology of [6, Sec. 3.6], which reduces to the stan- 
dard geometric distribution when po = 0. 

Conditionally on 6 and /i, the random variables X„, n > 1 
are independent; Xi, . . . , Xq-i and Xg^ ^e+i, ... are identi- 
cally distributed with common distributions Pq and P^, respec- 
tively. The probability measures Pq, Pi, . . . , Pjv/ always ad- 
mit densities with respect to some sigma-finite measure m on 
{E, £); for example, we can take m = Pg+Pi • • •+Pj\/. So, we 
fix m and denote the corresponding densities by /o , /i , . . . , Jm , 
respectively. 

Suppose now that we observe sequentially the random vari- 
ables Xn, n > 1. Their common pdf /o changes at stage 
to some other pdf /i G AA. Our objective is to detect the 
change time 9 as quickly as possible and to identify the change 
index /i, as accurately as possible. More precisely, given costs 
associated with detection delay, false alarm, and false identifi- 
cation of the change index, we seek a strategy that minimizes 
the expected total change detection and identification cost. 

Let F = {J-n)n>Q denote the natural filtration of the obser- 
vation process X, where 

= {0, ^} and Tn = cr(Xi, . . . , n > 1. 

A strategy 5 = (r, d) is a pair consisting of a stopping time t 
of the filtration F and a terminal decision rule d : ^ ^ 
measurable with respect to the history Tt = a{XnAT', n > 1) 
of observation process X through stage r. Applying a strategy 
6 = (r, d) consists of announcing at the end of stage t that the 
common pdf has changed from /q to fd at or before stage r. Let 

A = {(r, d) I r G F, and d G Tr is an X-valued r. v.} 

denote the collection of all such sequential decision strategies. 

For every strategy S =^ (r, d) G A, we define a Bayes risk 
function 

R{5) = cE[(r - 6)+] + E[ao<jl{^<e} + a^dl{e<r<oo}] (1) 

as the expected diagnosis cost: the sum of the expected detec- 
tion delay cost and the expected terminal decision cost upon 
alarm, where c > and > 0,z G {0} U A^, j G are 
known constants satisfying an = 0, z G (i.e., no cost for a 
correct terminal decision), and = max{a;, 0}. 

The problem is to find a sequential decision strategy 5 — 
(r, d) G A (if it exists) with the minimum Bayes risk 

R* 4 inf R{S). (2) 

5g A 



III. Posterior analysis and formulation as an 

OPTIMAL STOPPING PROBLEM 

In this section we show that the Bayes risk function in (1) 
can be written as the expected value of the running and termi- 
nal costs driven by a certain Markov process. We use this fact 
to recast the minimum Bayes risk in (2) as a Markov optimal 
stopping problem. 

Let us introduce the posterior probability processes 

n(°) ^ P{6i > n I Tn} and ^¥{9 <n,^^i\ Tn} 

for z G , n > 0. Having observed the first n observations, 
nl"' is the posterior probability that the change has not yet 
occurred at or before stage n, while n„ is the posterior joint 
probability that the change has occurred by stage n and that the 
hypothesis yu = i is correct. The connection of these posterior 
probabilities to the loss structure for our problem is established 
in the next proposition. 

Proposition 1. For every sequential decision strategy 5 G A, 
the Bayes risk function (1) can be expressed in terms of the pro- 
cess n = {n„ = {Un'' , . . . , ni*^'')}„>o as 



R{S) = E 



n=0 



M M 



i=0 



While our original formulation of the Bayes risk function (1) 
was in terms of the values of the unobservable random vari- 
ables 9 and /i. Proposition 1 gives us an equivalent version of 
the Bayes risk function in terms of the posterior distributions 
for 9 and /i. This is particularly effective in light of Proposition 
2, which we state with the aid of some additional notation that 
is referred to throughout the paper. Let 

^*'^{7r=(7ro,7ri,...,^M)G[0,l]^^+i 



1 



denote the standard M -dimensional probability simplex. De- 
fine the mappings Di : S^^ x E ^~> [0,1], i G M and 
D : S^' X E^ [0, 1] by 

, ^ f (1 -p)7ro.fo(2:), i = ] 
D^[n,x) = 



(tT; + TTo pvi)fi{x), i G M 



and D{tt, x) = X^fio -^il'"'; and the operator T on the col- 
lection of bounded functions / : 5*^^ h-^ R by 



(T/)W^/ ,n{dx)Din,x)f(^^ 
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for every tt G S^^ . 

Proposition 2. (a) The process II^") = {III"^ , JF„}„ 



(3) 



>o IS a su- 



permartingale, andW^Wn^ < (1 — pY^ for every n > 0. 

(b) The process II^*'' = {Iln\ J-n}n>o o. submartingale 
for every i G M. 



(c) The process 11 = 
process, and 

D,{nn,Xr 



{(n 



. . , Il''n^^)}n>o is a Markov 



D{Un,Xn+l) 



ie{0}UM, n>0, (4) 



with initial state IIq"'' ~ l—po and Hq' = p^Vi, i £ M.. More 



{i) 



over, for every bounded function f : S^'^ 

WE[/(n„+i)|n„] = (T/)(n„). 



. and n > 0, we 



Remark 3. Since 11 is uniformly bounded, the limit 
lim„_>oo n„ exists by the martingale convergence theorem. 
Moreover, lim„_too n^"'' = Q a.s. by Proposition 2(a) since 
(0,1). 

Now, let the functions h, hi, . . . , Hm from S^^ into IR+ be 
defined by 

AI 

h{Tr) = min hj{n) and hj{n) = Tij ajj , j G 

1=0 

respectively. Then, we note that for every 6 = (t, d) G A, we 
have 



R{t, d) = E 



> E 



M 
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.n=0 



where we define on the event {r < oo} the terminal deci- 
sion rule d to be any index satisfying h^{Ilr) = /i(ni-). In 
other words, an optimal terminal decision depends only upon 
the value of the 11 process at the stage in which we stop. Note 
also that the functions h and hi, ... , Hm are bounded on S^^ . 
Therefore, we have the following: 



Lemma 4. The minimum Bayes risk (2) reduces to the follow- 
ing optimal stopping of the Markov process II: 

R* = inf R{t, d) = inf R{t, d) 

{T.d)eA {rJjeA 



inf E 



^c(i-nW) + i{,<^}Mn,) 
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We simplify this formulation further by showing that it is 
enough to take the infimum over 

C:^{rGF|r<ooa.s. andEy- <cx)}, (5) 

where we define 

n-l 

-y„^^c(i-n["Vi(n,0, > o 

k=0 

as the minimum partial risk obtained by making the best termi- 
nal decision on {r = n}. Since /i(-) is bounded on S^^ , the 
process {Yn,J-n;n > 0} consists of integrable random vari- 
ables. So the expectation EYr exists for every t G F, and our 
problem becomes 



-R* = supEy^. 



(6) 



Observe that Er < oo for every t G C because oo > 

(l/c)Ey- > E(t-6I)+ > E{T-d) > Et-W > Et-(1/p). 
In fact, we have Ey^ > -oo <^ Ey^^ < oo <^ Er < oo for 
every r G F. Since sup^g^Ey^ > Eyo > ~h{Ilo) > -oo, it 



is enough to consider r G F such that Er < oo. Namely, (6) 
reduces to 



-i?* = supEy^. 

TEC 

IV. Solution via optimal stopping theory 



(7) 



In this section we derive an optimal solution for the prob- 
lem in (2) by building on the formulation of (7) via the tools of 
optimal stopping theory, which are detailed in [7] . 

A. The optimality equation. 

We begin by applying the method of truncation with a view 
of passing to the limit to arrive at the final result. Define for 
every pair of integers n, N satisfying < n < the sub- 
collections 

C„^{rVn|rGC} and = {t h N \t e C,,} 

of stopping times in C of (5) and the families of (truncated) 
optimal stopping problems 

(8) 



-Ki = sup Ey^ and - = sup Ey^ 



corresponding to (C„)„>o and {C^)Q<:n<N, respectively. Note 
that C = Co and R* = Vb- 

To investigate these optimal stopping problems, we intro- 
duce versions of the Snell envelope of {Yn)n>o (i-e., the small- 
est regular supermartingale dominating (y„)„>o) correspond- 
ing to (C„)„>o and (C^)o<«<Ar, respectively,'defined by 

7„ = ess supE[yT- I and 7,^ = ess supE[y^ | 

Then through the following series of lemmas we point out sev- 
eral useful properties of these Snell envelopes. Finally, we 
extend these results to an arbitrary initial state vector and es- 
tablish the optimality equation. Note that each of the ensuing 
(in)equalities between random variables are in the P-almost sure 
sense. 

First, these Snell envelopes provide the following alternative 
expressions for the optimal stopping problems introduced in (8) 
above. 

Lemma 5. For every N > and < n < N, we have 

-K = E7„flnfl'-V;^=E7,^. 

Second, we have the following backward-induction equa- 
tions. 

Lemma 6. We have 7„ = max{y„, E[7„-|_i \J-n]} for every 
n > 0. For every N > 1 and < n < N — I, we have 
In = and"f^ = max{y„, E[7^_|_i | Tn]}. 

We also have that these versions of the Snell envelopes coin- 
cide in the limit as iV ^ 00. That is. 



Lemma 7. For every n>Q, we have 7„ = lim 



N- 



' In ■ 



Next, recall from (3) and Proposition 2(c) the operator T and 
let us introduce the operator M on the collection of bounded 
functions / : i-^ M+ defined by 

(M/) W ^ min{;i(7r), c(l - ttq) + (T/)(7r)}, tt G . 



Observe that < Mf < h. That is, tt i-^ (M/)(7r) is a Proposition 15. For every N > 0, the function V^^ : 5^ 
nonnegative bounded function. Therefore, M^/ e 



IS continuous. 



well-defined. If / is nonnegative and bounded, then M"/ = 
M(M"- V) is defined for every n > 1, with M°f = f by defi- 
nition. Using operator M, we can express (7,^ )o<n<JV in terms 
of the process 11 as stated in the following lemma. 

Lemma 8. For every N > 0, and < n < N, we have 



N 



{0)^ 



1'. 

The next lemma shows how the optimal stopping problems 
can be rewritten in terms of the operator M. It also conveys 
the connection between the truncated optimal stopping prob- 
lems and the initial state Hq of the 11 process. 

Lemma 9. We have (a) Vq^ = (M'^/i)(no) /or every ^ > 0, 
and(b)Vo=^ lim (M^/i)(no). 

Observe that since Hq ^ J-q ~ {0, il}, we have P{no — 
tt} ~ 1 for some tt e S^^ . On the other hand, for every 
TT £ S^^ we can construct a probability space (fi, P^) host- 
ing a Markov process 11 with the same dynamics as in (4) and 
P^jllo = tt} = 1. Moreover, on such a probability space, 
the preceding results remain valid. So, let us denote by E^r the 
expectation with respect to Ptt and rewrite (8) as 



-Vn{n) = sup EttFt and 
Tec„ 



sup Ett^V 



for every tt e S . Then Lemma 9 implies that 



Vfl^TT) = (M'^ h){TT) and Vo{tt) = lim (M'^/i)(7r) (9) 

N^OQ 

for every tt e 5*^^. Taking limits as — > 00 of both sides 



m 



h){TT) ~ M(M /i)(7r) and applying the monotone 
convergence theorem on the right-hand side yields Vo{tt) = 
(MVb)(7r). Hence, we have shown the following result. 

Proposition 10 (Optimality equation). For every tt E S^^ , 
K)W = (MK))W =min{/iW,c(l-^o) + (Tyo)W}. (10) 

Remark 11. By solving Vo{tt) for any initial state tt € S^^, 
we capture the solution to the original problem since prop- 
erty (c) of Proposition 2 and (7) imply that R* = Vo(l ^ 
Po,PoVi, ■ ■ ■ ,PoVm)- 

B. Some properties of the value function. 

Now, we reveal some important properties of the value func- 
tion Vb( ) of (9). These results help us to establish an optimal 
solution for Vo(-), and hence an optimal solution for R* , in the 
next subsection. 

Lemma 12. If g : S^^ i-^ M /i a bounded concave function, 
then so is Tg. 

Proposition 13. The mappings tt V^{tt),N > and 
TT i-^ Vo{tt) are concave. 

Proposition 14. For every N > 1 and tt G we have 

\hr \\h\\\^l_ 
N' 

Since \\h\\ = sup^ggju \h{TT)\ < 00, limjv^oo i ^o^(^) 
Vq (tt) uniformly in tt G S^^. 



VoiTT) < < VoiTT) + ( ^ + ^ ) _. 



Corollary 16. The function Vq : S 



M 



1+ IS contmuous. 



Note that S^^ is a compact subset of R^^+^, so while conti- 
nuity of Vb(-) on the interior of follows from the concavity 
of Vb(-) by Proposition 12, Corollary 16 establishes continuity 
on all of S'^^, including its boundary. 

C. An optimal sequential decision strategy. 

Finally, we describe the optimal stopping region in im- 
plied by the value function Vo( ), and we present an optimal 
sequential decision strategy for our problem. Let us define for 
every > 0, 

T^^{7TeS'''\Vo''iTT)^hi7T)}, 
rj^' = n {tt G S^' I h{TT) ^ h^iTT)}, j G M, 

T ^ {tt eS''\ VoiTT) ^Htt)}, 

r(J) r n {tt G S''*^ | h{TT) = hj{TT)}, j g m. 

For each j G {0} U M, let ej G 5^^ denote the unit vec- 
tor consisting of zero in every component except for the jth 
component, which is equal to one. Note that eo, ■ • ■ , cm are 
the extreme points of the closed convex set , and any vec- 
tor TT ~ {ttq, . . . , ttm) G S^' can be expressed in terms of 

CM aSTT l^.^oTTjCj. 



eo, 



Theorem 17. For every j G Ai, (r^^)jv>o <^ decreasing se- 
quence of non-empty, closed, convex subsets of S^^ . Moreover, 

r(^) 3 {tt G S'^' I hj{TT) < min{/i(7r),c(l - ttq)}} 9 e^, 

00 A/ 00 

r = fi = y r(^), and r^-'") = fl m. 

N=l i = l N=l 

Furthermore, S^^ = To 3 Fi I) . . . I) F ^ {ei, . . . , cm}- 
Lemma 18. For every n > 0, we have 7„ = ^^^^^^^(1 — 

n("')-^o(n„). 

Theorem 19. Let a = inf{n > | n„ G F}. {a) The stopped 
process {7„Acrj J'n\ JT. > 0} is a martingale. 

(b) The random variable a is an optimal stopping time for 
Vq, and 

(C)E(T < QO. 

Therefore, the pair (ct, d* ) is an optimal sequential decision 
strategy for (2), where the optimal stopping rule a is given by 
Theorem 19, and, as in the proof of Lemma 4, the optimal ter- 
minal decision rule d* is given by 

d* — j on the event {a ~ n, Tin G F^^' } for every n >0. 

Accordingly, the set F is called the stopping region implied by 
Vo{-), and Theorem 17 reveals its basic structure. We demon- 
strate the use of these results in the numerical examples of Sec. 
V. 



Note that we can take a similar approach to prove that the 
stopping rules ctat = inf{n > | n„ G rjv-n},^ > are 
optimal for the truncated problems V(^{-), iV > in (9). Thus, 
for each > 0, the set F^r is called the stopping region for 
Vg'^(-): it is optimal to terminate the experiments in Fjv if ^ 
stages are left before truncation. 

V. Special cases and examples 

A. A. N. Shiryaev's sequential change detection problem. 

Set ooj = 1 for j e M and = for i,j e M, 
then the Bayes risk function (1) becomes R{5) = P{r < 
9} + cE[(r — 9)^]- This is the Bayes risk studied by Shiryaev 
[8, 9] to solve the sequential change detection problem. 

B. Sequential multi-hypothesis testing. 

Set pq = 1, then 6 = Q a.s. and thus the Bayes risk function 
(1) becomes R{5) ~ E[cr + a^dX{T<oo}]- This gives the se- 
quential multi-hypothesis testing problem studied by Wald and 
Wolfowitz [10] and Arrow, Blackwell, and Girshick [11]; see 
also [12]. 

C. Two alternatives after the change. 

In this subsection we consider the special case A/ = 2 in 
which we have only two possible change distributions, /i(-) 
and /2(-)- We describe a graphical representation of the stop- 
ping and continuation regions for an arbitrary instance of the 
special case M = 2. Then we use this representation to illus- 
trate geometrical properties of the optimal method (Sec. IV. C) 
via model instances for certain choices of the model parameters 
Po, p, vi, V2, fo{-), /i(-), /2(-)' aoi, ao2, ai2, 021, and c. 

Let the linear mapping L : M'^ be defined by 

L(7ro,7ri,7r2) = i^T^i + ^7i'2,7i'2)- Since ttq = 1 - tti - TTa 
for every tt — (ttq, tti, 772) G S"^ C M.'^, we can recover the 
preimage tt of any point L{tt) G L{S'^) C M^. For every point 
TT = (ttq , TTi , 772 ) G S'^ , the Coordinate tt^ is given by the Eu- 
clidean distance from the image point L(7r) to the edge of the 
image triangle L{S'^) that is opposite the image point L{ei), 
for each i = 0, 1, 2. For example, the distance from the image 
point L{tt) to the edge of the image triangle opposite the lower- 
left-hand corner L{1, 0, 0) = (0, 0) is the value of the preimage 
coordinate ttq. See Fig. ??. 

Therefore, we can work with the mappings i(F) and L{S'^ \ 
F) of the stopping region F and the continuation region 5^ \ F, 
respectively. Accordingly, we depict the decision region for 
each instance in this subsection using the two-dimensional rep- 
resentation as in the right-hand-side of Fig. ?? and we drop the 
L(-) notation when labeling various parts of each figure to em- 
phasize their source in S^. 

Each of the examples in this section have the following 
model parameters in common: 

f _('! 1 1 l^ f _/'J_ J3_ 2_ J_\ f _(J_ 2_ 3_ J_^ 
JO V4' 4' 4' 4/' V 10 ' 10 ' 10 ' 10/ ' •'2 V 10 ' 10 ' 10 ' 10 / • 

We vary the delay cost and false alarm/identification costs to 
illustrate certain geometrical properties of the continuation and 
stopping regions. See Figs. ??, ??, and ??. 



These figures have certain features in common. On each sub- 
figure there is a dashed line representing those states tt S S*^ at 
which /ii(7r) = h2{TT). Also, each subfigure shows a sample 
path of (nj,)^^g and the realizations of 9 and p for the sam- 
ple. The shaded area, including its solid boundary, represents 
the optimal stopping region, while the unshaded area represents 
the continuation region. 

Specifically, these figures show instances in which the M = 
2 convex subsets comprising the optimal stopping region are 
connected (Fig. ??) and instances in which they are not (Figs. 
?? and ??(a)). Fig. ??(b) shows an instance in which the con- 
tinuation region is disconnected. 

An implementation of the optimal strategy as described in 
Sec. IV.C is as follows: Initialize the statistic 11 = (n„)„>o 
by setting IIq = (1 ~ Po,PoI'i,Pqi^2) as in part (c) of Propo- 
sition 2. Use the dynamics of (4) to update the statistic n„ as 
each observation Xn is realized. Stop taking observations when 
the statistic n„ enters the stopping region F = F*^^^ U F*^^^ for 
the first time, possibly before the first observation is taken (i.e., 
n = 0). The optimal terminal decision is based upon whether 
the statistic n„ is in F^^^ or F^^^ upon stopping. Each of the 
sample paths in Figs. ??, ??, and ?? were generated via this al- 
gorithm. As Fig. ?? shows, the sets F^^^ and F^^^ can intersect 
on their boundaries and so it is possible to stop in their inter- 
section. In this case, either of the decisions d = 1 or d = 2 is 
optimal. 

We use value iteration of the optimality equation (10) over a 
fine discretization of to compute Vb( ) and generate the de- 
cision region for each subfigure. The resulting discretized deci- 
sion region is mapped into the plane via L. See [13, Ch. 3] for 
techniques of computing the value function via the optimality 
equation such as value iteration. 

D. Three alternatives after the change. 

In this subsection we consider the special case i\/ = 3 in 
which we have three possible change distributions, /i(-), /2(')' 
and /3(-). Here, the continuation and stopping regions are sub- 
sets of C K''. Similar to the two-alternatives case, we intro- 
duce the mapping of C M"* into via (ttq, tti, 7r2, tts) 

\\j2^^ + 2 V 2^2 + 2 V 2'^3, 2 \/ ^'^Z + 2 V 2^3,7^3^. 

Then we use this representation — actually a rotation of it — to 
illustrate in Fig. ?? an instance with the following model pa- 
rameters: 

PO-^, V^^V2 = Vz = \, 

JO V4'4'4'4/' Jl VIO' 10' 10' 10/ ' 

f. — (1. 1. JL ±.\ f — (3- 2- Jl. A.\ 

•'2 V 10 ' 10 ' 10 ' 10/ ' •'3 V 10 ' 10 ' 10 ' 10/ ' 

c=l, aoj = 40, fly = 20, i,j = 1,2,3. 

Fig. ?? can be interpreted in a manner similar to the fig- 
ures of the previous subsection. In this case, for every point 
TT = (ttq, TTi, 7r2, TTa) G , the coordinate tt; is given by the 
(Euclidean) distance from the image point L{'k) to the face of 



the image tetrahedron L{S^) that is opposite the image corner 
L{ei), for each i = 0, 1, 2, 3. 
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